Positron-atom scattering using pseudo-state energy shifts 
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A method to generate low-energy phase shifts for elastic scattering using bound-state calculations 
is applied to the problem of e + -Mg and e + -Zn scattering after an initial validation on the e + -Cu 
system. The energy shift between a small reference calculation and the largest possible configuration 
interaction calculation of the lowest energy pseudo-state is used to tune a semi-empirical optical 
potential. The potential was further fine-tuned by utilizing the energy of the second lowest pseudo- 
state. The s- and p-wave phase shifts for positron scattering from Mg and Zn are given from 
threshold to the first excitation threshold. The e + -Mg cross section has a prominent p-wave shape 
resonance at an energy of about 0.096 eV with a width of 0.106 eV. The peak cross section for 
e + -Mg scattering is about 4800 ag while Z c g achieves a value of 1310 at an energy of 0.109 eV. 
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One of the most technically demanding problems in 
quantum physics is the scattering problem, i.e. the pre- 
diction of the reaction probabilities when two objects col- 
lide The underlying difficulty lies in the unbounded 
nature of the wave function. This leads to a variety of 
computational and analytic complications that are ab- 
sent in bound state calculations, e.g. the Schwartz sin- 
gularities that occur in the Kohn variational method for 
scattering 0,(3]. 

One approach to solve scattering problems is to use 
methods that have been used for bound state calculations 
[E 0, HI]. There are many examples of such approaches, 
one of the most popular being the i?-matrix methods 
that use the solutions of the Schrodinger equation in a 
finite sized cavity to determine the behavior of the wave 
function in the interaction region [J [H, 0, HI • The total 
wave function is constructed by splicing the inner wave 
function onto the asymptotic wave function. 

This article had its origin in a particular scattering 
problem, namely the determination of the near thresh- 
old phase shifts for positron scattering from the di-valent 
group II and IIB atoms. The dimension of the secular 
equations for bound state calculations on such systems 
are very large, for example a configuration interaction 
(CI) calculation of the e + Ca 2 P° state resulted in equa- 
tions of dimension 874,448 These dimensions are 
much larger than those that would occur in a CI calcula- 
tion of the 2 P° ground state Ca~ . The exceptionally large 
dimensionalities occur because the valence electrons tend 
to localize around the positron, thus giving a very slowly 
convergent partial wave expansion of the wave function 
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PH El G3, El Q d. The ability to routinely solve 
the secular equations associated with the CI basis us- 
ing iterative sparse matrix techniques [HI is one reason 
why CI calculations for positronic atoms (and of course 
for molecular systems) have been able to generate useful 
results. 

Trying to generate scattering solutions for such sys- 
tems would be problematic for a number of reasons. For 
example, application of the CI-Kohn approach [f| to de- 
termine the phase shifts for positron scattering from any 
group II or IIB atom would result in linear equations 
with dimensions between 400,000 to 1,000,000. These are 
simply too large to be solved by direct methods. Itera- 
tive methods for large linear systems do exist, but there 
are no robust methods that absolutely guarantee conver- 
gence [IT}- The development of an efficient linear solver 
for the class of problems that arise from a basis set treat- 
ment of quantum scattering would likely involve a good 
deal of initial experimentation and effort. Similarly, the 
widely used i?-matrix method with fixed boundary con- 
ditions (3] requires the generation of all the eigenvectors 
and eigenvalues of the Hamiltonian, which is not feasible 
when the matrix dimensions exceed 100,000. 

Very recently, a method was developed to extract phase 
shifts from the positive energies of a pair of CI calcula- 
tions [l8j]. In that work, the energy shifts of a positive 
energy pseudo-state were used to tune a semi-empirical 
optical potential which was then used to predict the close 
to threshold phase shifts. This concept is refined in the 
present work and it is shown that the reliability of the 
potential can be enhanced by tuning the potential to the 
energies of the two lowest states. Next, the s- and p-wave 
phase shifts for positron scattering from Mg and Zn are 
computed from threshold to the opening of the lowest 
energy excitation channel. The Mg and Zn atoms are 
interesting for positron scattering experiments since Mg 



2 



has been recently shown to possess a prominent p-wave 
shape resonance [Hj]. The Zn system is also interesting 
since the existence of a e + Zn bound state of 2 S C symme- 
try [l9| will manifest itself in a differential cross section 
that is largest at backward angles [20| . 



I. MODEL INDEPENDENT METHOD FOR 
GENERATING PHASE SHIFTS 

A. The box variational method 

The idea behind the current method lies closest to the 
box variational method [2l|, [22|, [23| which is exploited 
in quantum Monte Carlo (QMC) calculations of scatter- 
ing [3, [H, [H, H3, El- In the box variational method, 
one extracts the phase shift by comparing the zero point 
energy of a finite size cavity to the energy of the sys- 
tem wave function in the same cavity. In its simplest 
incarnation for s-wave scattering, one diagonalizes the 
Hamiltonian, in natural units where Ti = m = e = 1, 
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V(r), 



(1) 



in a cavity of radius R. The wave function obeys the 
boundary conditions ^(0) = *£(R) = 0. The positive 
energy states $ n (r), with energy E n , can be regarded as 
the small r-part of the exact scattering wave function, 
ty n (f), with that same energy. The exact wave function 
can be written as ^ n {r) — sin(fc„r + S) for r > R where 
5 is the phase shift and the wavenumber k n — y/2E n . At 
the boundary, one has sin(k n R + S) = 0, giving 



S = mr — k n R, 



(2) 



(this expression assumes there are no bound states). 

For systems with non-zero angular momentum, the 
asymptotic sin-wave is replaced by the asymptotic form 
ip(r) ~ jt{kr) + tan(<5j»)rifc(fcr) where ji(kr) and ne(kr) 
are spherical Bessel functions of the 1st and 2nd kind. 
The condition ip(R) — gives the following expression 
for the phase shift, 



tan(^(fc„)) = 



je(k n R) 

m(k n R) 



(3) 



B. Phase shifts using pseudo-state energy shifts. 

The box variational method has two advantages, (a) 
it is very simple to apply and (b) the -B-spline basis sets 
currently in use in many atomic structure applications 
easily satisfy the necessary boundary conditions. How- 
ever, there are other basis sets in use that do not satisfy 
the ^(0) = ^(R) = boundary conditions. 

Consider the usage of a set of general L 2 functions, u — 
{4>i}. These functions have a finite radial extent and thus 
the basis can be regarded as defining a soft-sided cavity. 



A simple procedure is used here to estimate the radius of 
the resulting soft-sided box. Denoting E[u,0] n to be the 
nth energy eigenstate resulting from a diagonalization of 
the V = potential in the basis u, then the effective 
radius of the soft box is given by 
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X, 
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^2E[u,Q] n k[u,0} n 



(4) 



where Xi n is the nth zero of the spherical Bessel function, 
k{x). 

The potential V — V(r) is then diagonalized in the 
same basis to give E[u,V] n , and hence k[u, V] n — 
y / 2E[u, V] n . The phase shift is then extracted using 



je(k[u,V] n R) 
ni(k[u, V] n R) 



(5) 



Figure[T]shows the radial probability density, ^(r)] 2 = 
I c a r Xa\ 2 , of the four lowest energy I = wave func- 
tions computed by the diagonalization of a V — poten- 
tial in a basis of 30 Laguerre Type Orbitals (LTOs) with 
the scale parameter A = 1. The general definition for the 
LTOs were 

Xa{r) = N a r e exp(-A Q r)4 2 ^ 1 (2A a r) , (6) 
where the normalization constant is 



l (2X a ){n a -i-l)\ 



D! 



(7) 



The function L^+|l 1 (2A a r) is an associated Laguerre 
polynomial that can be defined in terms of a confluent 
hyper-geometric function [29| . The probability densities 
go to zero for r > 60. The wiggles in the probability 
densities are not a numerical artifact, rather they are a 
manifestation of the slow convergence of the L 2 basis to 
the exact continuum wave function 301. 
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FIG. 1: The radial probability densities, ^(r)! 2 , for the four 
lowest s-wave pseudo-states resulting from a diagonalization 
of the V = potential in a basis of 30 LTOs. The pseudo- 
states have all been normalized to unity. 



3 



Table U gives the energies, and the effective radius, of 
the soft box for this basis as given by eq. (pf]). The four 
states have an effective radius of about 61.5 which is con- 
sistent with Figure [1] It is reasonable to conclude that 
eq. d3|) gives an estimate of the range of a pseudo-state 
that is sufficiently accurate to be useful. 

This LTO basis was also used to diagonalize the 
Woods-Saxon potential, 



V(r) 



V 



1 + exp 
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with the choice Vq = 0.97, W = 1 and a = 0.05. The 
energies and phase shifts derived from eq. (O are listed 
in Table [IJ The phase shifts obtained by numerically in- 
tegrating the Schrodinger equation for the Woods-Saxon 
potential are also listed in Table U and are exact to all 
quoted digits. The two sets of phase shifts agree with 
each other to an accuracy of about 2%. 



FIG. 2: The radial probability densities, |^(r)| 2 , for the four 
lowest p-wave pseudo-states resulting from a diagonalization 
of the V = potential in a basis of 30 LTOs. The pseudo- 
states have all been normalized to unity. 



II. MODEL DEPENDENT METHOD FOR 
GENERATING PHASE SHIFTS 



TABLE I: Parameters derived from the diagonalization of the 
free-wave and Woods-Saxon potential in a basis of 30 LTOs 
with A = 1. The pseudo-state energies for V = is denoted 
Eq while the Woods-Saxon energies are denoted i? ws . The 
radius of the soft-box is denoted Rq, while the phase shift 
from eq. (0) is 5. The phase shift obtained by integrating the 
Woods-Saxon potential numerically is <5 C xact- 



n 


E 


Ro 


-Ews 


6 


f^exact 








£ = 






1 


0.001286 


61.95 


0.001170 


0.1451 


0.1470 


2 


0.005170 


61.79 


0.004705 


0.2897 


0.2852 


3 


0.011734 


61.52 


0.010736 


0.4098 


0.4092 


4 


0.021116 


61.15 


0.019362 


0.5333 


0.5148 








1= 1 






1 


0.002474 


63.87 


0.002450 


0.0212 


0.0211 


2 


0.007361 


63.67 


0.007147 


0.1114 


0.1103 


3 


0.014810 


63.36 


0.013964 


0.3132 


0.3073 


4 


0.024969 


62.94 


0.022844 


0.6087 


0.5903 



This procedure has also been validated for p-wave scat- 
tering. Figure [2] shows the result of diagonalizing the 
V = potential for p-wave scattering in a basis of 30 
LTOs with 1 = 1 and A = 1.0. Once again the range of 
the pseudo-state solutions are roughly the same. Table 
H] lists the effective box radius for each pseudo-state as 
derived from eq. ((4|). 

The Woods-Saxon potential with the choice Vq = 
0.173, Wo =4.0 and a = 1.0 was then diagonalized in 
this basis. The phase shifts obtained from eq. (O are 
tabulated in Table U along with phase shifts generated 
by a numerical solution of the Schrodinger equation. The 
two sets of phase shifts agree to within 3%. 



It has long been known that the positive energy 
pseudo-states resulting from the diagonalization of the 
Hamiltonian in an L 2 basis often give a reasonable ap- 
proximation to the exact scatte ring wave function over a 
finite range [H HI H [H, [H [36| . Methods which ex- 
ploit this result are sometimes called stabilization meth- 
ods. While convergence of the pseudo-state to the con- 
tinuum wave function is relatively slow, matrix elements 
formed by the pseudo-state often have reasonable conver- 
gence properties [30|. In effect, while the point- wise prop- 
erties of the wave function can be inaccurate, the conver- 
gence in the mean of the wave function over a suitable 
range can be quite good. This raises the possibility that 
the expectation values of positive energy pseudo-states 
can be used to define a semi-empirical optical potential 
to describe low-energy scattering. 

Our method proceeds as follows. The initial calcu- 
lation uses a reference basis of square integrable single 
particle orbitals, {^(r)}, designed to give a good repre- 
sentation of the wave function in a bounded interaction 
region. The Hamiltonian, Hq for a free particle with 
V = is diagonalized, yielding the wave function 

$o = 5> s: <Mr) , (9) 

i 

and the energy expectation, 

E bee = ($ |#o|$o>. (10) 

The wave function of the target atom is constructed 
from a linear combination of configurations {aj;(X)}, of 
the same symmetry as the ground state (X is the collec- 
tive set of target coordinates). So one can write 

r> gs (x) = ^ Cl ^(x) , (ii) 
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while 

E gs = (r2gs|-fftarget|^gs)- (12) 

The reference energy, E , is determined by diagonalizing 
the Hamiltonian in the product basis, f2 gs (X)0j(r), to 
give 

E Q = (* |fftargct|*0> , (13) 

where 

* = 5><Mr)ft gs (X) ■ (14) 

t 

The basis sets {4>i(r)} and {wi(X)} are then aug- 
mented by a large number of additional functions {xi( r )} 
and {^(X)} to represent the correlations between the 
projectile and the target constituents. None of these 
additional functions have the same symmetries as those 
used in {^(r)} and {cj, (X)}. This augmented trial func- 
tion can be written as 

*i = S c^(X)^(r) + ]T rf i>J -Vi(X) Xj (r) . (15) 

The trial wave function, v&i, is used to diagonalize H cxact 
giving an energy 

E X = (^lliJcxactl*!)- (16) 

Next, the basis {(f>i(r)} is diagonalized in a parameter- 
ized potential designed to describe the most important 
features of the interaction between the projectile and the 
target. This potential can be written formally as 

V opt (r)=V dh (r)+V po i(r) . (17) 

The potential Vdi r is the direct interaction between the 
target and projectile. This can be approximated by the 
direct interaction between the projectile and the tar- 
get Hartree-Fock (HF) ground state wave function, IIhf, 
which can be slightly different from f2 gs . The polariza- 
tion potential V po i(r) is semi-empirical in nature with the 
asymptotic form 

^poi(r) ~ -|J , (18) 

where ay is the static dipole polarizability of the target. 
In previous work, [l8l |20| . r3TL \3^ a simple one-parameter 
form 

F P i(r) = -^(l-exp(-r 6 /p 6 )), (19) 

has usually been adopted for V po i(r). It is thought that 
this functional form has the incorrect shape at interme- 
diate values of r, e.g. r ss 5 ao (The reasons why we orig- 
inally become suspicious about the reliability of eq. (fl"9|) 
are not discussed here. But the results obtained later will 
clearly show the limitations of this type of cutoff polar- 
ization potential) . The present work will also use a more 



complicated expression for V po i(r), with an additional ad- 
justable parameter, Aq to give an improved description 
of the potential between target and atom. This form was 

V p2 (r) = -gL(l-exp(-r 6 /p 6 )) 

- ||(l_exp(-r 8 /p 8 )). (20) 

While the second term has the functional form of a 
quadrupole polarization potential, it should be regarded 
as primarily empirical in nature. This functional form 
was chosen as a screened quadrupole type potential be- 
cause it was computationally convenient. 

The energy expectation value of the ground state, or 
lowest energy pseudo-state 

E opt = <$opt|ffopt|$opt}, (21) 

is adjusted by tuning the parameters in V p i until E op t — 
Ei — E gs . 

Determination of the V P 2 required additional informa- 
tion since there are two parameters, p and Aq that need 
to be fixed. In this case, the optical potential is tuned to 
two energy levels rather than one. This does increase the 
overall time of the calculation since it is necessary to ex- 
tract the lowest two eigenvalues from the CI calculation. 

Once the optical potential has been fixed, it is a sim- 
ple matter to generate the exact continuum solution of 
the Schrodinger equation for the Hamiltonian given by 
eq. mU). 



A. Positron annihilation 

Besides obtaining the phase shifts in the low-energy 
region, it is also possible to determine the annihilation 
parameter, Z e g [lCj, |39j, l40( • The fundamental idea is to 
compare exact and model potential calculations of Z e g, 
and so fix the enhancement factor, G [H, |4l|, [42j. En- 
hancement factors were first introduced in the calculation 
of the annihilation rate of positrons in condensed matter 
systems [43|, 0, ■ They incorporate the tendency for 
attractive electron-positron correlations to increase the 
electron density at the position of the positron. 

It has been shown that model potential calculations 
of s-wave positron scattering from hydrogen and helium 
that were tuned to give the correct phase-shift at a refer- 
ence energy also reproduced the low-energy behavior of 
Z e g(k) up to a multiplying constant (i.e. G) [HI]. The an- 
nihilation parameter for the model potential wave func- 
tion follows the model of Mitroy and Ivanov [38| , and is 
written as 

Z cfi = J d 3 r(G v p v (r)+G cPc (rj^ |$ opt (r)| 2 , (22) 

where p c ( r ) an d Pv( r ) are the electron densities associ- 
ated with the core and valence electrons of the target 
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atom, and <J> opt (r) is the positron scattering function ob- 

tained in the tuned model potential. The notation Z c g 
is used to denote the annihilation parameter for the £th 
partial wave. 

For the core orbitals, G c is set to 2.5 due to reasons 
outlined in Ref. 38] . The valence enhancement factor G v 
is computed by the simple ratio 



where is the annihilation rate of the positron with 
the valence orbitals as given by the CI calculation and 
pmodei j g |_ ne va i ence annihilation rate predicted by the 
model potential calculation with G = 1. 

III. THE FIXED CORE POTENTIALS 

All calculations on the e + Cu, e + Mg, and e + Zn sys- 
tems used a fixed core Hamiltonian. The details of the 
corepotentials have been discussed previously 0, HE 
EE S3) El]> but a short description is worthwhile. The 
model Hamiltonian is initially based on a HF wave func- 
tion for the neutral atom ground state. One and two- 
body semi-empirical polarization potentials are added to 
the potential field of the HF core and the parameters of 
the core-polarization potentials defined by reference to 
the spectra of Cu, Mg+ and Zn+ [H HE EE S3 • 

The effective Hamiltonian for the systems with 2 va- 
lence electrons (r^ and r 2 ) and a positron (r ) was 

1 2 1 

H = -2 V »^2 V '~ l ' di,(ro) + l ' cpl(ro) 

2—1 

2 2 x 

+ YAKlir(ri) + V cxc (r t ) + V cp i{Ti)) - V — 
i=l i=l n ° 

1 2 

+ V r C p2(ri,r2)+2^y c p2(ri,ro) . (24) 

^ i=l 

The direct potential (Vdir) represents the interaction with 
the HF electron core. The direct part of the core po- 
tential is attractive for electrons and repulsive for the 
positron. The exchange potential (V^c) between the va- 
lence electrons and the HF core was computed without 
approximation. 

The one-body core polarization potentials (V cp i) are 
semi-empirical in nature. They have the functional form 

K P i(r)=-£^fPlM<M- (25) 

The factor ad is the static dipole polarizability of the 
core and gf{r) is a cut-off function designed to make the 
polarization potential finite at the origin. The same cut- 
off function has been adopted for both the positron and 
electrons. In this work, gj(r) was defined to be 

. 9 f(r) = l-exp(-r 6 /^) , (26) 



where pi is an adjustable parameter. The two-body po- 
larization potential (V^ p2 ) is defined as 

V r cp2 (r l , Vj) = -g-g- (i-i • r 3 )g cp 2(n)g cp 2(r 3 ) . (27) 

i 3 

where g C p2{r) is chosen to have a cut-off parameter, p cp 2, 
obtained by averaging the pi. The core dipole polariz- 
abilities were set to 0.4814 for Mg [HE1, 5.36 for 
Cu dzj], and 2.294 for Zn The cutoff parameters 
for Mg were p = 1.1795 a , pi = 1.302 a , p 2 = 1-442 
a , P3 = 1.52 a , pe>4 = 1-361 a , and p cp2 = 1.361 
a . The cutoff parameters for Cu were po — 1-9883 a , 
p x = 2.03 a , p-2 = 1-83 a , p 3 = 1.80 a , p^> 4 = 1.91 a , 
and p cp 2 — 1.91 ciq. The cutoff parameters for Zn were 
po = 1.63 a , pi = 1.80 a , p 2 = 2.30 a , p 3 = 1-60 a , 
pe>i = 1-83 ao, and p C p2 = 1-83 clq. This model has been 
used to describe many of the features of neutral Be, Mg, 
Ca and Sr to quite high accuracy [IE EE EH • 



IV. VERIFICATION FOR e+-Cu SCATTERING 

Previously, a validation of the method was performed 
for s-wave e + -H scattering fl8| . In the present work the 
method is further verified by computing the low-energy 
phase shifts and annihilation parameters for s-wave and 
p-wave e + -Cu scattering. The model copper atom used 
here has a dipole polarizability of 41.65 [29J and, there- 
fore, provides a more stringent test of the procedure used 
to tune the shape of the polarization potential than the 
previous test upon the e + -H system (where ad = 4.5 q,q). 

The explicit CI calculation on the e + -Cu ground state 
and the CI-Kohn calculations of e + -Cu scattering closely 
follow those previously reported [5j, |47|, [48| • Briefly, the 
wave function expansion consists of a large number of 
single particle orbitals and includes terms with I > 10. 
The single particle orbitals are usually represented as La- 
guerre type orbitals (LTOs). 

The e + Cu ground state calculation included orbitals 
up to L < 16 with a minimum of 18 electron LTOs and 
18 positron LTOs per £. The CI reference wave function, 
^0) consisted of the copper atom ground state multiplied 
by a positron basis of 30 £ = LTOs. The orbital ba- 
sis was slightly reduced for the calculation of the lowest 
energy 2 P° pseudo-state. In this case, the calculation in- 
cluded terms up to L — 14 with a minimum of 18 electron 
LTOs and 18 positron LTOs per £. The CI reference wave 
function, ^0) i n this case consisted of the copper atom 
ground state multiplied by a positron basis of 33 £ = 1 
orbitals. 

One difficulty present in all CI calculations of positron- 
atom interactions is the slow convergence of the energy 
with L [3, El) HO]- The con verg ence pattern of the 
atomic CI expansion [H, [U, HE I53l. l54l HE HH , suggests 
the use of an asymptotic analysis that utilizes the result 
that successive increments, AEl = (E)l — (E)l-i, can 



TABLE II: Expectations values obtained from CI calculations 
of the bound 2 S e states and the 2 P° pseudo-states for the 
e+Cu, e + Mg and e + Zn systems. The binding energy e (in 
Hartree) is negative for bound states and positive for pseudo- 
states. The mean positron radius (r p ) is in units of ao. The 
core ((r) c ) and valence ((r„)) annihilation rates are given in 
units of 10 9 s _1 . All of the values given in this table are the 
results of extrapolating L — > oo. 



System 


Symmetry 


£ 


(r P ) 


<r} c 




e+Cu 


2ge 


— 0.005124 


9.037 


0.0322 


0.5035 


e+Cu 


2po 


0.001860 


35.23 


0.000413 


0.0186 


e+Mg 


2ge 


-0.01704 


6.930 


0.0109 


1.004 


e+Mg 


2po 


0.003989 


13.87 


0.00110 


0.3729 


e+Zn 


2 S e 


-0.003794 


9.726 


0.0244 


0.4269 


e+Zn 


2po 


0.006885 


20.24 


0.000609 


0.0190 



eq. (123)) was evaluated for the ground state. The details 
of the model potential parameters are summarized in Ta- 
ble ETTl 
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be written as an inverse power series, viz. 

Be Ce 



AE r 



A, 



+ 



(28) 



The L — > oo limits have been determined by fitting sets 
of (E)l values to asymptotic series with either 1, 2 or 3 
terms. The factors, Ae, Be and Ce for the 3-term ex- 
pansion are determined at a particular L from 4 succes- 
sive energies ((E) L -3, (E)l-2, (E)l-i and (E) L ). The 
series is summed to oo once the linear factors have been 
determined. 

Some expectation values of the e+Cu 2 S e ground state 
and the lowest energy 2 P° pseudo-state in the L — > oo 
limits are given in Table [TTJ It should be noted that the 
leading term of the inverse power series for the annihi- 
lation rate, T, is A r /(L + 1/2) 2 [H, HI]. There is some 
uncertainty in the extrapolation procedure and we esti- 
mate uncertainties of about 1% for the energy and 5% 
for the annihilation rate. However, this does not im- 
pact the present verification exercise. The extrapolation 
procedures were applied consistently to both the CI cal- 
culations used to define the model potentials (and en- 
hancement factors) , as well as the independent CI-Kohn 
scattering calculations [j| used to validate the model po- 
tential calculations. Note that the errors in the extrap- 
olated results introduced by the use of a finite basis set 
have a tendency to fortuitously cancel out [4&| . 

The trial function, was then used to diagonalize 
the model potential, eq. ([2TJ)1 with two different polariza- 
tion potentials. In the first instance, eq. (Til?]) was used 
and the parameter p varied until the energy matched that 
of the CI calculation. This potential will be referred to 
as the V p i potential. In the second instance, the param- 
eters, Aq and p of eq. (f2T))) were both varied until both 
the energy of the ground state and lowest energy pseudo- 
state were the same as the CI calculations. This potential 
will be termed the V P 2 potential. The enhancement fac- 
tor, G„, was determined after the model potentials were 
finalized. In the case of the V P 2 potential the ratio in 



FIG. 3: Comparison of the two different parameterizations of 
the e+-Cu 2 S e polarization potential of eq. (|19[) and eq. (120[) 
against that of the asymptotic form of eq. (|18|l . 

Figure [3] shows a comparison of the V p \ and V P 2 po- 
larization potentials for the 2 S e symmetry. The V p i po- 
tential is always smaller in magnitude than the ad/(2r 4 ) 
asymptotic form. The V P 2 potential bulges below the 
ad/(2r 4 ) asymptotic form and is stronger than a pure 
dipole potential in the outer valence region of the atom. 
This is entirely reasonable. The slow convergence of the 
single-center expansion occurs as a result of the localiza- 
tion of the valence electrons in the vicinity of positron 
[H ED, ES El Q El- This, in turn, enhances the 
strength of the polarization potential in the outer valence 
region. 

The superiority of the V P 2 potential in describing the 
the 2 S e bound state is apparent from Tables ITT1 and UTTl 
The V p 2 calculation overestimates the core annihilation 
rate by 26% while the V p i potential overestimates this 
parameter by more than 120%. Additionally, the V P 2 
potential gives a better estimate of the mean positron 
radius, (r p ). The value of 8.822 ao is about 2% smaller 
than the CI value of 9.037 ao while V p i gave (r p ) — 8.46 
a (~ 6% smaller). 

Accurate phase shifts for the full e+-Cu scattering 
Hamiltonian were obtained from CI-Kohn variational cal- 
culations [f| of the e+-Cu system using exactly the same 
short-range orbital basis sets as used in the CI calcu- 
lation. The only difference between the CI-Kohn and 
regular CI basis sets is the addition of two continuum 
basis functions [5j|- The phase shifts for the V p i and V P 2 
potentials were obtained by integrating the Schrodinger 
equation. 

The V P 2 scattering length estimate of 12.8 ao is within 
2% of the CI-Kohn estimate of the scattering length, 
namely 13.05 ao- The V^i scattering length of 12.4 ao is 
too small by 5%. Figure S] shows the comparison between 
the model potential s-wave phase shift and the CI-Kohn 
phase shift for k S [0,0.2] a^ 1 . The V p \ model slightly 
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TABLE III: Definitions of the V p i and V P 2 model potentials used to describe s-wave and p-wave scattering of e + -Cu, e + -Mg 
and e + -Zn. The s-wave potentials were tuned to the properties of the 2 S e ground state and the lowest energy pseudo-state, 
while the p-wave potentials were tuned to the two lowest energy p-wave pseudo-states. The binding energy e (in Hartree) is 
negative for bound states and positive for pseudo-states. The mean positron radius and scattering length, j4 SC at, are in units 
of ao- The core and valence annihilation rates are given in units of 10 9 s" 1 . 



Atom 


Potential 


L 




Aq 


P 


G v 


e 


(r v ) 


<r c > 


(T„) 




Cu 


V P i 





41.65 


0.0 


2.7434 


18.94 


-0.005124 


8.46 


0.0730 


0.5036 


12.4 


Cu 







41.65 


480.0 


3.6248 


26.35 


-0.005124 


8.822 


0.04088 


0.5034 


12.8 


Cu 


v pl 


1 


41.65 


0.0 


2.1231 


20.01 


0.0057801 


35.18 


0.00275 


0.0186 




Cu 


V p2 


1 


41.65 


360.0 


3.0829 


36.70 


0.0057801 


35.22 


0.000868 


0.0186 




Mg 







71.35 


0.0 


2.9927 


13.12 


-0.017072 


6.21 


0.0243 


1.004 


6.09 


Mg 


V p2 





71.35 


2280.0 


4.4794 


24.74 


-0.017072 


6.982 


0.00738 


1.004 


7.23 


Mg 


Vpi 


1 


71.35 


0.0 


2.5626 


12.35 


0.003989 


12.90 


0.00654 


0.3729 




Mg 


v p2 


1 


71.35 


1250.0 


3.8406 


28.15 


0.003989 


13.80 


0.00115 


0.3729 




Zn 


Vpi 





41.25 


0.0 


2.6579 


9.91 


-0.003794 


9.34 


0.0412 


0.4269 


14.3 


Zn 


V p2 





41.25 


430.0 


3.5344 


14.35 


-0.003794 


9.71 


0.0219 


0.4269 


14.7 


Zn 


Vpi 


1 


41.25 


0.0 


2.1604 


10.45 


0.006885 


20.21 


0.00177 


0.0190 




Zn 


V p2 


1 


41.25 


252.0 


3.0117 


17.45 


0.006885 


20.26 


0.000770 


0.0190 





overestimates the CI-Kohn phase shifts (modulo w) over 
the entire range. The V P 2 fit to the CI-Kohn phase shifts 
is clearly superior. 

Besides obtaining phase shifts, this procedure was used 
to determine the valence annihilation parameter which is 
shown in Figures [6] and [Jj The V P 2 enhancement factor 
of G v — 26.35, gives an s-wave annihilation parameter, 

zjS , that is within 5% of the explicit CI-Kohn calcula- 
tion over the entire energy range. Somewhat surprisingly, 
the from the Vpi model is almost the same as that 
from the V" P 2 model. 




FIG. 4: The phase shift for e -Cu scattering in the s-wave as 
a function of k (in units of a^ 1 ). The lines show the results of 
the present calculation using the tuned Vpi and V P 2 potentials 
while the squares show the phase shifts of the explicit CI-Kohn 
calculation. 

The phase shift for p-wave scattering is shown in Figure 
[5j The Vpi potential overestimates the CI-Kohn phase 
shift as the energy increases and there is a 15% discrep- 



ancy at k = 0.20 a Q . The Vpi potential also tends to 
underestimate the phase shift for k < 0.10 a^ 1 , although 
this is difficult to see from the Figure. The V P 2 potential 
reproduces the CI-Kohn phase shifts very well and the 
agreement is perfect within the resolution of the graph. 
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FIG. 5: The phase shift for e -Cu scattering in the p-wave as 
a function of k (in units of a^ 1 ). The lines show the results of 
the present calculation using the tuned V p i and V^,2 potentials 
while the squares show the phase shifts of the explicit CI-Kohn 
calculation. 

This pattern is repeated in Figure [7] where is plot- 
ted as a function of k. The Vpi potential tends to overes- 
timate the CI-Kohn values at the higher momenta with 
the discrepancy at k = 0.20 a^ 1 being 15%. However, 
the Vp2 potential does an excellent job of reproducing 

the CI-Kohn Z^g over the entire momentum range. The 
Vp2 ^off i s t°° large at the higher momentum, but the 
difference is only 2% at k = 0.20 Oq . 
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FIG. 6: The annihilation parameter, for e + -Cu s-wave 
scattering as a function of k (in units of a^ 1 ). The two curves 
were calculated with the Vpi and V P 2 potentials. The discrete 
points are taken from the explicit CI-Kohn calculations. 
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7: The annihilation parameter for p-wave scattering, 
for e + -Cu scattering as a function of k (in units of a^ 1 ). 
The different curves were calculated with the V p i and V P 2 
potentials. The discrete points are taken from the explicit 
CI-Kohn calculations. 



The final confirmation of the improved quality of the 
V P 2 potential comes from the comparisons of the pseudo- 
state expectation values in Tables |TT] and The (r p ) 
given by V V 2 is closer to the CI value than the V p i value. 
Furthermore, the V P 2 potential is better than the V p \ 
potential at reproducing the CI core annihilation rate of 
0.0322 x 10 9 s" 1 (this value assumes G c = 1). 



V. POSITRON SCATTERING FROM 
MAGNESIUM 

A. The CI calculations 

Although many of the specifics of the calculations upon 
e + Mg have been reported previously [H, further de- 
tails concerning the wave function construction are given 



here. The trial wave function adopted for the CI calcula- 
tions consists of a linear combination of states which are 
anti-symmetric in the interchange of the two electrons, 



\V;LS) a = J2 c i\®i-> LS )a 



(29) 



Each anti-symmetrized state is constructed as a linear 
combination of coupled but not anti-symmetrized states. 
Two electrons (particles 1 and 2) are coupled first to each 
other, then the positron (particle 0) is coupled to form 
a state with net angular and spin angular momentum, L 
and S. The anti-symmetric states are written as 



[ab]L I S I pLS) a 



y/2(l + S ab ) 
+ (-l) eML < +s >\[a 2 b 1 ]L I S m )) , 

(30) 

where the subscript by each orbital denotes the electron 
occupying that particular orbital. 

The e + Mg CI basis was constructed by letting the two 
electrons and the positron form all the possible configu- 
ration with a total angular momentum of L, with the two 
electrons in a spin-singlet state, subject to three selection 
rules, 



max(£ ,4,^2) < L , 
min(^i,£ 2 ) < L int , 



= +1 or - 1 . 



(31) 
(32) 
(33) 



In these rules lo, l\ and £2 are respectively the orbital 
angular momenta of the positron and the two electrons. 
The even [odd] parity states require (— l)( £ o+ £ i+^) = 

[-!]• 

The Hamiltonian for the e + Mg 2 S° state was diagonal- 
ized in a CI basis including orbitals up to £ = 12. There 
were a minimum of 15 radial basis functions for each I. 
There were 19 i = positron orbitals. The largest 2 S e 
calculation was performed with L = 12 and L- mt = 4. 
The Li n t parameter does not have to be large since it is 
mainly concerned with describing the more quickly con- 
verging electron-electron correlations The CI ba- 
sis for the e + Mg 2 P° symmetry included orbitals up to 
£ = 14. There were a minimum of 14 radial basis func- 
tions for each £. There were 20 £ — 1 positron orbitals. 
The largest 2 P° calculation was performed with L = 14 
and Lint = 3. 

A summary of e + Mg expectation values taken to the 
L — > 00 limit are given in Table [TTJ The binding en- 
ergy £ for each symmetry is calculated with respect 
to the energy of the Mg ground state using the basis 
for that symmetry. The overall binding energy of the 
2 S e ground state was —0.017072 Hartree, with the first 
pseudo-state at 0.002503 Hartree. The energies of the 
two lowest pseudo-states of 2 P° symmetry were 0.003989 
and 0.012012 Hartree respectively. 
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B. Model potential calculations 
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FIG. 9: The elastic scattering cross section for e + -Mg scatter- 
ing as calculated with the V P 2 potential in the energy region 
FIG. 8: The s, p and d-wave phase shifts for e+-Mg scattering below the Ps-formation threshold at k « 0.249 a 1 . The solid 



as a function of k (in units of a^ 1 ). The solid lines were 
computed with V P 2 while the dashed lines were computed with 
V p i. The horizontal line shows 5 — tt/2. 



The p parameter of the V p i potential was tuned to re- 
produce the energies of the lowest 2 S e and 2 P° states. 
The values of p and the expectation values of the lowest 
state of each symmetry are given in Table ITlTl The values 
of Aq and p for the V P 2 potential were tuned to the low- 
est two energies. Examination of the expectation values 
of Tables |TT] and IIIII reveals that the V P 2 potential again 
does a better job at reproducing the CI expectation val- 
ues. The Vpi potential underestimates the mean positron 
radius by 10% and further overestimates the core anni- 
hilation rate by a factor of 2. The V P 2 potential gives a 
value of (r p ) that is too large by 1%. The V P 2 underesti- 
mation of T c is about 30%. While the V P 2 model potential 
may not be perfect, it does a better job of describing the 
radial distribution of the positron density than the V p i 
potential. 

The situation for the 2 P° pseudo-state is similar to that 
for the 2 S e state. The V p \ potential underestimates the 
mean positron radius by 10% and overestimates the core 
annihilation rate by a factor of 6. The V P 2 potential on 
the other hand gives an (r p ) within 1% of the CI value 
and overestimates the core annihilation rate by only 5%. 

The s- and p-wave phase shifts are plotted in Figure 
13 The 10% difference between the two model potential 
scattering lengths manifests itself in the slightly differ- 
ent s-wave phase shifts. The difference between the V p \ 
and V P 2 potentials is larger for the p-wave phase shift, 
although both predict a resonance at k w 0.09 O.Q 1 . The 
(i-wave phase shift plotted in Figure[5]was computed with 
the V P 2 p-wave potential. The £ > 2 phase shifts used in 
the computation of the total cross section also used the 
V P 2 p-w&ve potential. 

Figure [5] shows the elastic cross section for e + -Mg scat- 
tering below the Ps formation threshold (at k w 0.25 a^ 1 ) 
as computed with the V P 2 potentials. The p-wave reso- 



line shows the total cross section while the dashed curve shows 
the £ — 1 partial cross section. 
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The existence of the resonance in the phase shift also 
leads to a resonance in the annihilation parameter, Z e s 
[l~8| . The curve in Figure [TU] was computed with an en- 
hancement factor of G c = 2.5 for core annihilation, an 
s-wave valence enhancement factor of G v — 24.7 and a 
p-wave enhancement factor of G v — 28.2. The total Z c s 
shown in Figure |TD] is almost completely dominated by 
the contribution from the p-wave and the value of Z c g at 
the resonance peak was 1310. 

The resonance parameters were determined by per- 
forming a fit to the function 



S = S + a(E — er) + tan 



F 



2(e R - E) 



bE 2 . (34) 
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This gave a value of er = 0.00351 Hartree for the reso- 
nance position and a width of T — 0.00390 Hartree. 



C. Reliability of the resonance prediction 

The possible sources of error in the e + -Mg calcula- 
tions are (a) the reliability of the underlying model po- 
tential for the CI calculation, (b) the extent to which 
the CI calculations have converged and (c) the ability of 
the scattering model potential to reliably reproduce the 
scattering parameters. 

The cutoff parameter for the positron part of the core 
polarization potential (refer to eq. (|2"5")1 ) is chosen to be 
the same as the cutoff for the electron. This is likely 
to underestimate the strength of the positron interaction 
since there is a good deal of evidence for closed shell sys- 
tems that suggests the positronic part of the polarization 
potential is stronger than the electronic part [501 ] . How- 
ever, the impact of this is likely to be small since the 
core polarizability of 0.4814 aj] is more than 100 times 
smaller than the neutral atom polarizability of 71.35 a^. 
Any correction would tend to shift the resonance to a 
lower energy and increase the height of the maximum in 
the cross section. 

The CI calculations of the 2 S e ground state are believed 
to be converged to about 2% in the energy. An inde- 
pendent calculation of the e + Mg ground state has been 
done with the fixed core stochastic variational method 
(FCSVM) The FCSVM Hamiltonian is very sim- 

ilar to the fixed core Hamiltonian used for the present 
calculation and the current best FCSVM estimate of the 
binding energy is 0.017117 Hartree. However, it has also 
proved possible to make an estimate of the variational 
limit of the FCSVM calculation. This estimate is be- 
tween 0.01735 and 0.01740 Hartree [H] which is about 
2% more tightly bound than the CI calculation. The cal- 
culation of the 2 P° state is expected to have an accuracy 
similar to that of the ground state. 

The existence and position of the resonance is indepen- 
dent of the exact form of V po i. Besides the calculations 
reported here, alternate calculations with some other pa- 
rameterizations were reported earlier [l8j]. All of these 
calculations gave a resonance almost at the same po- 
sition and magnitude. The reason for this lies in the 
accident that the energy of the 2 P° pseudo-state, at 
k = 0.0893 a,Q , lies close to the center of the resonance. 
At this energy, the determination of the phase shift will 
be largely model independent since the stabilization con- 
cept ensures that the L 2 wave function is a reasonable 
approximation to the actual continuum wave function. 
The phase shifts of the two potentials at the energy of 
the pseudo-state were 8\ = 1.157 and 5\ = 1.153 radi- 
ans, for Vpi and V P 2 respectively. Additional plots of 
the phase shifts obtained with other functional forms for 
V^oi have tended to have a common intersection point 
near k w 0.089 a^ 1 . 

Finally, the simple potential independent approach of 



eq. ([5]) has been applied to determine the phase shift at 
the pseudo-state energy. The energy of the positron p- 
wave LTO basis in the V — potential was 0.007572 
Hartree. The radius of the box giving this energy is 
-Rbox = 36.5 ao- Evaluating eq. ([5]) at k — 0.0893 a^" 1 
gives S = 1.16 radian. 



VI. POSITRON SCATTERING FROM ZINC 

Positron binding to zinc has been known with some de- 
gree of certainty since 199 9 fl9l.[2 0| following some earlier, 
less conclusive, work [59l 1 6 0L 1 6 ll] . The neutral zinc atom 
has an ionization potential of 0.34523 Hartree |62| and a 
polarizability of 38.8±0.8 ag [63]. The present model po- 
tential for the Zn 2+ core predicts an ionization potential 
of 0.33519 Hartree and a polarizability of 41.25 tig [201 ]. 

The present CI calculations upon the e + Zn ground 
state used the same core potential as the earlier CI cal- 
culations [20], but the size of the basis has been en- 
larged. The maximum number of partial waves has been 
increased to L — 12, the number of LTOs per £ has been 
increased to 16, and finally L- mt was increased from 3 to 4. 
The overall dimension of the CI calculation has increased 
by an order of magnitude. The summary of e + Zn expec- 
tation values for the series of calculations with increasing 
L are given in Table IIVI The energy of the Zn ground 
state with respect to the Zn 2+ core for the electron basis 
was -0.99549251 Hartree. 

An examination of Table IIVI reveals that the present 
extrapolated binding energy of 0.0037944 Hartree is 
within 2% of the previously obtained binding energy. To 
a certain extent, this high level of agreement is fortu- 
itous. The method used to extrapolate the energy incre- 
ment to the L — > oo limit in Ref. (20| had an inherent 
tendency to overestimate the binding energy. However, 
this compensated for the tendency of a finite dimension 
LTO basis to increasingly underestimate the energy in- 
crement as L — ► oo [48]. The lowest positive energy 2 S e 
pseudo-state had an energy of 0.0024706 Hartree above 
the Zn ground state. 

The Hamiltonian for the e + Zn 2 P° state was diagonal- 
ized in a CI basis including orbitals up to £ — 10. The 
two electrons were in a spin-singlet state, with a mini- 
mum of 16 radial basis functions for each £. There were 
20 £ = 1 positron orbitals. The largest calculation was 
performed with L = 10 and L m t = 3. The energy of 
the lowest energy 2 P° pseudo-state was 0.006885 Hartree 
above threshold. Other expectation values for this state 
are listed in Table HTl The second lowest 2 P° pseudo-state 
was located at 0.019055 Hartree. 

The parameters for the V p i and V P 2 potentials, tuned 
to the CI data in Table [TT1 are listed in Table [TTTJ A ca- 
sual glance at the entries in these two Tables reveals that 
the V P 2 potential again does better at reproducing the 
properties of the 2 S e physical state and the 2 P° pseudo- 
state. 

The s- and p-wave phase shifts are plotted in Figure 
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TABLE IV: Results of the CI calculations for e + Zn atoms for a given L. The E column gives the energy with respect to the 
doubly ionized frozen core and e is the binding energy with respect to the lowest energy dissociation channel at E — —0.99549251 
Hartree. The radial expectation values (in ao) of the electron and positron are listed in the (r e ) and (r p ) columns. The (r„) 
and (r c ) columns give the valence and core annihilation rates (in 10 9 sec -1 ). The results in the row labeled 10* are taken from 
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FIG. 11: The s- and p-wave phase shifts for e + -Zn scattering 
as a function of k (in units of ). The solid lines were 
computed with V P 2 while the dashed lines were computed with 
Vpi. The horizontal dashed line shows 5 — tt/2. 



1111 The small e + Zn binding energy leads to a large value 
for the scattering length, namely 14.7 ± 0.1 ao. The en- 
ergy region below k < 0.144aQ 1 would be interesting for 
an experimental investigation as the s-wave phase shifts 
are in a different quadrant from all the other phase shifts. 
Consequently, the differential cross section will be larger 
at backward angles than at forward angles. This peaking 
of the near threshold differential cross section at back- 
ward angles is a signature of the existence of the e + Zn 
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FIG. 12: The elastic scattering cross section for e -Zn scat- 
tering as calculated with the V P 2 potential in the energy region 
below the Ps-formation threshold at k ~ 0.436 a^ 1 . The solid 
line shows the total cross section while the two dashed curves 
show the £ = and £ = 1 partial cross sections. 



bound state (201 ] - 

The low-energy elastic cross section below the Ps for- 
mation threshold as shown in Figure [T^] was computed 
with the V P 2 potentials. The phase shifts for £ > 1 are 
taken from the 1 = 1 model potential. The large value of 
the cross section at E = is characteristic of a potential 
supporting a weak bound state. The quickly rising p- 
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FIG. 13: The annihilation parameter Z c g for e + -Zn scatter- 
ing as a function of k (in units of a^ 1 ). The different curves 
show, Z^g , Ztf{ and ^^ otal ) as calculated with the V P 2 po- 
tential 

wave phase shifts leads to a shoulder in the cross section 
near k « 0.2 Oq . The p-wave causes a more pronounced 
structure in Z c g which is easily noticeable in figure IT51 as 
the bump at k » 0.2 a^ 1 . 

The experimental observation of the e + Zn p-wave reso- 
nance precursor in the cross section would be complicated 
by the large s-wave cross section which tends to obscure 
the feature in the total elastic cross section. The reso- 
nance structure would be most visible in a measurement 
of Z c g or in a differential cross section. 



VII. COMPARISONS WITH PREVIOUS WORK 

The present calculations are not the only calculations 
of the e + Mg and e + Zn scattering systems. However, 
the other calculations were of a much more speculative 
nature [13, H H M, HE HI M, Isl- For example, 
the many body perturbation theory-based calculation of 
Gribakin and King predicted that the 2 P° symmetry of 
e + Mg had a bound state (66[. None of the other calcula- 
tions on the e + -Mg system gave a cross section with the 
prominent 2 P° shape resonance. 

While some previous model potential calculations were 
based on reasonable estimates of the e + Mg binding en- 
ergy [53, HH , the uncertainties in defining the functional 
form of the polarization potential detracted from the reli- 
ability of the p-wave phase shift. The present calculations 
are more definitive, and the main source of uncertainty 
is in the definition of the underlying core polarization 
potential in the CI calculations. 



VIII. SUMMARY 

A new technique has been used to determine the phase 
shifts for low-energy positron-atom scattering from mag- 



nesium and zinc. The phase shifts are determined by tun- 
ing an optical potential to the energy of a bound state or 
a positive energy state. The tuning of an optical potential 
to features such as bound state energies and resonance 
positions is well known. Tuning an optical potential to a 
pseudo-state energy shift is novel [1 81 ] . One improvement 
over our previous implementation is the use of a second 
energy to fine-tune the shape of the optical potential. 
Another possible improvement requiring further research 
would be to use other expectation values (e.g. the mean 
positron radius) to further refine the shape of the optical 
potential. The use of experimental information to tune 
optical potentials is known (e.g. the role of the deuteron 
radius in tuning the n-p potential). 

There are two different concepts that can be regarded 
as providing motivation for the present approach. The 
first is the stabilization concept, namely, a positive en- 
ergy pseudo-state will provide a reasonable approxima- 
tion to the scatte ring state with that energy over a finite 
radial range [3TI . |32| . The alternate motivation comes 
from the box variational method, namely that the energy 
shift of the wave function in a hard-sided box is used to 
estimate the phase shift [U [22|, [23| . Diagonalization of 
the Hamiltonian in a finite dimension LTO basis can be 
regarded as equivalent to diagonalizing the Hamiltonian 
in a soft-sided box. 



T 1 — I I I I I I 1 1 1 — I I I I I I 1 1 1 — I I I I I 




0.01 0.1 1 10 

Energy (units of eV) 

FIG. 14: The cross sections (in 10 3 a§) for e + -Mg scattering 
versus energy (in eV). The present elastic scattering cross 
sections as calculated with the V P 2 potential are shown as 
the solid line. The total cross section measurements from the 
Detroit group arc shown as the discrete points with their error 
bars [63I . The inset shows the same data, plotted on a linear 
energy scale. 

The most significant result of the present investigation 
is the prediction of a close to threshold 2 P° shape res- 
onance for elastic scattering from magnesium. There is 
no experimental evidence for the existence of shape reso- 
nances in positron-atom or positron-molecule scattering 
[o0| . The present prediction has the virtue of being read- 
ily amenable to experimental verification. Indeed, the 
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Detroit group has measured the total cross section for 
positron-magnesium scattering down to an energy of 2.0 
eV . As can be seen from Figure [HJ the lowest en- 
ergy for which their measurements were done is just too 
high to detect the resonance. An earlier experiment mea- 
sured down to an energy of 1.0 eV [zj|, but these results 
are not shown in the figure as they are similar to those 
in Ref. (69j whilst having larger reported errors. Their 
most recent e + -Mg measurements went down to 0.12 eV, 
however, they only reported the positronium formation 
cross section (7l| . 

With some reflection on the differences between 
positron-atom and electron-atom interaction potentials 
it is not surprising that the magnesium atom supports a 
shape resonance. It has been noticed that the positrons 
are more strongly attracted to closed (sub)shell atoms 
than are electrons (5(| (this result is based on results for 
systems with 2 S e symmetry). Since there is a low-energy 
2 P° shape resonance in e~-Mg scattering [72|, [73[, one 
could reasonably infer that the e + -Mg system would also 
have a 2 P° shape resonance or, alternatively, support a 
bound state. 

The e + -Zn cross section has a broad feature in the p- 
wave at about 0.6 eV that could be interpreted as a res- 
onance or a precursor to a resonance. It should be noted 
that a similar feature occurs in e~-Zn scattering 74| at 
roughly the same energy. Figure [15] compares the p-wave 
phase shift from the B-spline i?-matrix (BSRM) calcu- 
lation of e~-Zn scattering [74j with the present e + -Zn 
phase shift. While the BSRM is probably not converged 
with respect to the enlargement of the channel space, the 
low-energy elastic cross section does a reasonable job at 
reproducing the electron transmission experiment of Bur- 
row et al [751 ]. The similarity between the electron and 
positron p- wave phase shifts for k < 0.10 a^ 1 is expected 
since the low-energy phase shifts will be dominated by 
the long range polarization potentials. It is interesting 
to speculate upon whether the polarization potential will 
lead to e + -Zn phase shifts that are larger than the e~-Zn 
phase shifts for the 2 P° symmetry, as well as the 2 S e sym- 
metry. The comparison depicted in Figure [T5l shows that 
the e~-Zn 2 P° phase shift is larger than the e + -Zn phase 
shift for k > 0.14 . This is possibly due to the elec- 
tron seeing an attractive static potential as it penetrates 
the centrifugal barrier while the positron experiences a 
repulsive potential. However, it would be best to test 
this conjecture using models for the Zn target which are 
exactly the same. 

The existence of the e + -Mg resonance and the struc- 
ture in e + -Zn suggest that other group II atoms might 
support a 2 P° shape resonance. The dipole polarizability 
of beryllium is only a bit smaller than that of Zn, so some 
sort of structure in the 2 P° partial wave is expected. The 
cadmium atom, on the other hand has a larger polariz- 



ability than zinc, so a more pronounced resonance should 
be expected. 

The actual polarization potentials used here represent 
a departure from those used in some previous calculations 




k (units of a ) 



FIG. 15: The p-wave phase shifts as a function of k (in 



units of a 



o 



for 



-Zn scattering as given by the present 



calculation and the BSRM calculation [74 



of positron-atom interactions |20|, [3J, [3J, [JJ, [77j . All of 
these previous works use polarization potentials with a 
cut-off function that leads to a potential that is always 
smaller in magnitude than that of the ad/(2r 4 ) asymp- 
totic potential. The present polarization potentials have 
bulges in the outer valence region that are larger in mag- 
nitude than the asymptotic potential in that region. 

A possible area of application of the current approach 
would be to positron-molecule scattering. However, 
this would require improvements in the technology of 
positron-molecule CI calculations. The best calculations 
so far carried out [ll, [7^, HO] do not treat the electron- 
positron dynamics nearly was well as the present CI cal- 
culations on atoms. 
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